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Abstract 

In this letter we report results from inelastic neutron scattering exper- 
iments on powder samples of vanadyl pyrophosphate, (VO)2P207. We see 
evidence for three magnetic excitations, at 3.5 meV, 6.0 meV and 14 meV. 
The intensity of the 3.5 meV mode is strong near Q = 0.8 A -1 , consistent 
with the one-magnon gap mode reported previously and predicted by the 
spin-ladder model. The magnetic scattering at 14 meV may be due to the 
top of the one-magnon band or the two-magnon continuum predicted in the 
ladder model. The 6.0 meV mode also peaks in intensity at Q = 0.8A -1 . 
This mode has not been reported previously and was not anticipated by 
existing theoretical treatments of the spin ladder. 



I. INTRODUCTION 



The antiferromagnetic insulator (VO)2P207 [l| (abbreviated VOPO) is widely con- 
sidered to be the prototypical realization of a Heisenberg spin ladder. Heisenberg spin 
ladders are interesting theoretically M as intermediaries between ID Heisenberg antifer- 
romagnets and 2D systems. Theoretical studies of this model have concluded that S=l/2 
spin ladders have a gap if they have an even number of chains, but are gapless for an odd 
number of chains. This is reminiscent of the Haldane result that half-integer-spin chains 
are gapless but integer-spin chains have energy gaps. Studies of hole-doped 2-chain lad- 
ders in the t-J model find that strong d-wave hole pairing takes place on equal-strength 
ladders, so analogues of the high-Tc superconductors might be found in hole-doped spin 
ladders. Additional examples of spin ladders are known, including a Sr-Cu-0 series [|J, a 
La-Cu-0 series |4| and (^(CsH^^^C^ ||. The subject of spin ladders was reviewed 
recently by Dagotto and Rice 0. 

The magnetic properties of VOPO have been interpreted in terms of the 2-chain Hamil- 
tonian 

chains rungs 

h = J,, J $t ■ s, + j± E ^ ■ Sj . (i) 

<ij> <ij> 

This model gives an accurate fit to the VOPO magnetic susceptibility J7] with the param- 
eters J || = Jj_ = 7.8 meV. Theoretical studies of the ladder model predict an energy gap of 
Egap = 0.5037J H, hence the susceptibility results imply E sap = 3.9 meV. A pulsed neu- 
tron scattering experiment on VOPO powder || observed a gap of 3.7(2) meV, consistent 
with this expectation. 

Although these aspects of VOPO are consistent with the Hamiltonian (1), many unan- 
swered questions remain concerning the dispersion relations, spectral weight, and detailed 
band structure of the magnetic excitations. We have carried out triple-axis inelastic neu- 
tron scattering studies of VOPO to investigate these issues in more detail. In addition to 
confirming the existence of the gapped spin wave mode observed previously, our results 



2 



reveal a new and unexpected magnetic excitation. 

The spectrum of low-lying spin waves in the 2-chain ladder model has been discussed 
at length in the literature P,|7|JT0| . The one-magnon band can be approximately described 
by the function J/J 

r n 1/2 

u(k) = [uj{0) 2 cos 2 (£;/2) + cj(tt) 2 sin 2 (A;/2) + c 2 , sin 2 (A;)J (2) 

where k is the one-dimensional wavevector in units in which the intrachain ion separation 
is set to unity. This dispersion relation for the one-magnon band u(k) is shown in Fig.l 
together with the lower edges of the two-magnon and three-magnon continua. (These 
are the lowest-energy sums uj{k\) + uo(k2) and uj(k{) + u(k2) + uj(k^) with the constraint 
J2i ki = k.) The value J=7.0 meV was chosen to match the energy gap observed in the 
present work, as discussed below, and our dispersion relation parameters for Jj_ = Jy = J 
are w(0) = 1.72J, u>(tt) = 0.5037J and Cq = 1.38J; these values were motivated by finite 
lattice results. 

The lowest-lying one-magnon excitation occurs at the one- dimensional antiferromag- 
netic point k = ir, with an energy of 00(71) = E gap . Exact calculations of the dynamical 
structure factor S(k, uj) on a 2 x 8 lattice 0,(11]] indicate that the spectral weight of the 



one-magnon band peaks strongly near k = tt, so this should be the most prominent feature 
seen in neutron scattering. 

The two-magnon states span a continuum of levels, beginning at 2E gap at k = 0. 
As seen in Fig.l, the lower edge of the two-magnon continuum crosses the one-magnon 
band at an intermediate k w 0.37T, and for smaller k these two-magnon states are the 
lowest-lying excitations. Similarly there is a three-magnon continuum, also shown in 
Fig.l, which contains the first k = tt excitation above E 9ap in the ladder model; this 
three-magnon continuum begins at 3E gap . 

Numerical calculations of the structure factor S(k, u) |TT| indicate that the two- 
magnon states have moderate spectral weight near k = tt, at E « 2.2J (w 15 meV in 
VOPO), and the one-magnon band may also be observable near k = 0, at E w 1.7J 



(« 12 meV in VOPO). Both excitations have S(k, uj) ~ 0.1 £(& = 7r, = E 9ap ) of the gap 
mode. The S(k,u) of the three-magnon continuum at k = ir, E w 1.5J (~ 11 meV) is 
about a factor of five weaker than these higher excitations. 

II. EXPERIMENTAL DETAILS 

To prepare (VO)2P207, stochiometric amounts of V2O5 (44.1 grams, 99.99%) and 
NH4H2PO4 (55.9 grams, 99.9%) powder were mixed and loaded into a platinum crucible. 
The mixed powders were heated to 500°C at l°C/minute and held at 500°C for 6 hours. 
This initial heating allowed the ammonia and water to escape from the crucible during 
the decomposition of NH4H2PO4 to P2O5. The prereacted mixture was heated further to 
a temperature of 1100°C and the resulting liquid was allowed to homogenize for several 
hours. The platinum crucible was removed from the furnace at 1100°C and the liquid 
was "cast" into a cold, shallow platinum crucible measuring 10 x 4 x 2 cm 3 . The rapid 
cooling of the V205-P20s liquid resulted in the formation of a homogenous glass with a 
dark green to black color. The above steps were all performed in air. 

The shallow platinum crucible containing the V2CVP2O5 glass was then placed in a 
tube furnace in which the oxygen partial pressure could be controlled. A mixture of 0.1% 
oxygen and 99.9% argon was passed through the furnace at 40 cc/minute, and this rate 
and oxygen partial pressure were maintained during the entire crystal growth procedure. 
An independent oxygen sensor was used to monitor the oxygen partial pressure in the 
furnace. The V205-P 2 05 glass was heated to 1050°C, maintained at that temperature for 
6 h, and then cooled at l°C/hr to 480°C, after which the furnace was turned off. The 
resulting material consisted of several hundred mm sized single crystals of (VO)2?207 
embedded in a polycrystalline matrix of (VO)2P20y and a small amount (~ 5-15%) of an 
additional phase, tentatively identified as V(POs)3 by neutron powder diffraction [ p~2]j . 

Inelastic neutron scattering measurements were carried out using the HB1A and HB3 
triple-axis spectrometers at the High Flux Isotope Reactor at Oak Ridge National Lab- 
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oratory (T^j. The sample consisted of 18 grams of finely ground VOPO powder sealed 
in a 0.5" diameter thin-walled Al cylinder under a helium atmosphere. For the initial 
experiment on the HB1A spectrometer the sample was mounted in a standard He closed- 
cycle refrigerator with a temperature range of 10K to 300K. A double-crystal pyrolytic 
graphite (PG) (002) monochromator (M) was used to provide an incident neutron beam 
with a fixed energy of 14.7 meV. The incident beam intensity was monitored using a fission 
counter, and contamination of the beam by higher-order Bragg reflections was removed 
using the standard PG filter method. Neutrons scattered by the sample (S) were reflected 
by a PG(002) analyzer (A) into a 3 He detector (D). Effective beam collimations pre-M, 
M-S, S-A and A-D of 40'-40'-40'-70' or 40'-20'-20'-70' were used. 

A second experiment utilizing the HB3 spectrometer employed PG(002) monochro- 
mator and analyzer crystals with a fixed scattered neutron energy of 14.7 meV. The 
collimation was 40 / -40'-40'-120', and the scattered beam was passed through a PG fil- 
ter. For this experiment the sample was mounted in a pumped He cryostat capable of 
reaching 1.5K. 

VOPO has a slightly monoclinic unit cell which has lattice parameters a = 7.728 A, 
b = 16.589 A, c = 9.580 A, and (3 = 89.98° at T=296K ]TJ. The ladder is composed 
of S= 1/2 V 4+ ions in -(VO)-(VO)- chains along the a axis and V ZqZ V runes oriented 
along the b axis. The average V-V distance along the chain is approximately equal to 
a/2. Therefore in VOPO the one- dimensional antiferromagnetic points k = mr (with n 
odd) occur in sheets with Q a = nQ n , where Q n = 0.8134 A -1 . 

In scattering from powder the experimental momentum transfer to the sample Q = \Q\ 
does not correspond to the one-dimensional momentum k of the excitations observed. 
Instead k is the projection of Q on the ladder axis, and in the powder one averages over 
all orientations. Therefore for a given Q one can excite all ladder modes with the observed 
energy transfer and k < Q. In particular, one expects the gap mode, which has k = tt, to 
appear first at Q = Q n as Q is increased, and to persist to higher Q due to the powder 
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average. The scattering intensity observed at (Q,u) involves the density of states as well 
as S(k, u). 

III. RESULTS 

Figure 2 shows the results of constant-Q scans at Q = Q^, taken at HB1A. The data 
have been scaled to a fixed monitor value and have been corrected for the constant-Ej 
kinematical factor kf cot 9a- The data at 10. 8K clearly show two peaks in the scattering. 
The peak intensities are drastically reduced at higher temperatures, as illustrated by the 
data at 41K, which provides strong evidence that they are both magnetic in origin. The 
lower peak at 3.48(3) meV is interpreted as the one-magnon gap mode previously reported 
at 3.7(2) meV by Eccleston et al. 0. Their slightly higher energy may arise from their 
broader Q acceptance. The 3.5 meV peak observed in the present work has a FWHM 
significantly greater than the instrumental resolution limit of 0.7 meV. 

The mode at E = 6.0 meV has not been reported previously. It appears narrower in 
energy than the 3.5 meV mode, with a FWHM only slightly larger than the instrumental 
resolution; this suggests that the band is rather flat, possibly corresponding to a relatively 
localized mode. Comparison with the theoretical magnon bands in Fig.l shows that this 
excitation is not predicted by the simple Heisenberg spin-ladder model (1). 

Figure 3 shows the Q dependence of the intensities of the 3.5 meV and 6.0 meV modes. 
Both modes show very similar behavior, with an onset and peak at Q = Q^. This is pre- 
cisely what one would expect for the lowest-energy one-magnon gap mode. The fact that 
the 6.0 meV mode peaks at a value of Q characteristic of VOPO provides strong evidence 
that it actually is a VOPO excitation, and does not arise from a contaminant phase in the 
sample. Both modes also show a slight increase in intensity in the vicinity of Q = 2Q 7T , 
and a gradual falloff in intensity at higher Q. The scattering intensity from a powder sam- 
ple is not easily interpreted beyond the first Brillouin zone, but the diminishing intensity 
with increasing momentum transfer is expected because of the reduction of the magnetic 
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form factor at large Q. For the V 4+ ion in VOPO, \f ma9 {3Q«)\ 2 « 0.5|/ mos (Q 7r )| 2 . 

Constant Q scans at Q — 3Q n and Q = 5Q n at a temperature of 1.6K are presented in 
the lower panel of Fig. 4. The scans were carried out at HB3 with fixed scattered neutron 
energy, which allowed measurements over a wide range of energy transfers. The scattering 
at both wavevectors exhibits a broad maximum near 13 meV. The origin of this broad 
scattering is unknown, but since magnetic scattering at Q = 5Q n is expected to be weak, 
it is probably not magnetic in origin. Unfortunately, kinematic restrictions did not allow 
the same energy scan at lower Q, particularly at Q = Q n where the magnetic features 
should be most prominent. Inspection of the lower panel of Fig.4 shows that the 3Q- K scan 
has features which were not observed in the 5Q n scan. In repetitions (not shown) of the 
3Q tt scan at T=45K and T=100K the extra features are not clearly visible, and only the 
broad maximum remains. The Q and T dependence of the scattering strongly suggests 
that the additional features seen at 3Q- K arise from magnetic excitations. To investigate 
this possibility further the Q = 5Q n intensity was subtracted from the Q = 3Q n intensity 
on a point-by-point basis; the result is shown in the upper panel of Fig.4. The solid line is 
a fit of this difference spectrum to a sum of three Gaussians and a constant background, 
revealing features near 3.6 meV and 6.5 meV, approximately consistent with the magnetic 
excitations observed at Q n . The fit also shows a higher-energy peak at 14.3 meV. Eccleston 
et al. || also saw indications of enhanced scattering near 15 meV at Q = 0.8 A -1 , although 
the details were unclear. 

Since 14 meV is close to the predicted onset of the k = it two-magnon continuum as 
well as the maximum energy of the one-magnon band (see Fig.l) in the ladder model, 
it is tempting to associate the magnetic scattering with either or both of these features. 
However, a detailed comparison of experiment with the excitation spectrum predicted by 
(1) is problematical, since it did not anticipate the 6 meV mode. 



7 



IV. DISCUSSION 



In this experiment we have found clear evidence for two low-lying magnetic excitations 
in VOPO, at energies of approximately 3.5 meV and 6.0 meV, and some additional mag- 
netic scattering near 14 meV. The low-lying excitations peak strongly at Q = 0.8 A -1 , 
near the k = it point in VOPO, which suggests that these are k = ir modes of the VOPO 
spin ladder. Although the 3.5 meV mode and the 14 meV scattering are similar to ex- 
pectations in the spin-ladder model, the presence of two low-lying modes at k = it is 
unexpected, and appears to argue against the validity of the simple spin-ladder model (1) 
for VOPO. 

Of course one possibility is that the theoretical predictions we associate with the spin 
ladder, such as the spectrum and continua shown in Fig.l, are inaccurate. Although it 
appears highly unlikely that two low-lying excitations arise at k = ir in the model, these 
conclusions are based on Lanczos studies of relatively small clusters (up to 2 x 12), and 
careful theoretical investigations of larger spin ladders would be useful. 

Assuming that the theoretical interpretation of the model is accurate, we conclude 
that there are magnetic excitations in VOPO which are not consistent with the Heisenberg 
spin-ladder Hamiltonian (1). In view of this apparent disagreement we should recall the 
assumptions which led to (1) and determine whether they are justified. The Heisenberg 
model assumes that each magnetic ion can be treated as a localized spin with S= 1/2, 
interacting with neighboring ions through an isotropic interaction. Whether or not this is 
true in practice depends on the ground state and low-lying excited states of the magnetic 
ion. 

An isolated V 4+ ion in an octahedral field, as is approximately realized by the O 
octahedra in VOPO, has a t 2g triplet ground state [f[5| ; this triplet is split by the spin- 



orbit interaction and by departures from octahedral symmetry, including the formation of 
(VO) units along the chains fL4f . For isolated V 4+ ions in AI2O3 the single-ion excited t2 9 



levels are known to lie just 3.5 meV and 6.6 meV above the ground state |16[], comparable 
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to the energy scales of the magnetic modes observed in VOPO. 

This raises the possibility that one or more of the magnetic modes observed in VOPO 
could be crystal field levels, which might explain the peak at 6 meV. There is experimental 
evidence against this possibility, however; a single-ion excitation should not display a 
strong dependence of the scattering intensity on Q (other than the magnetic form factor). 
Instead we have observed a strong peak at Q n in the constant-E scan (Fig. 3), which 
suggests that the 6 meV excitation is part of a band of delocalized ladder excitations. 
Of course this band could involve excited t2 g orbital levels, in which case VOPO should 
be described by a multiband Hubbard model with t 2g V 4+ orbitals interacting through 
intermediate O levels. 

Although low-lying t 2g V 4+ orbitals may indeed be present in VOPO, magnetic sus- 
ceptibility measurements argue against their presence at these low energy scales. The 
observed high temperature limit of the susceptibility (known to about 350K [|IJ) is con- 
sistent with that expected for paramagnetic spins with S=l/2 and g=2. The presence 
of additional t 2g degrees of freedom would give a much larger susceptibility as the tem- 
perature approached the orbital gap. Thus there is no indication of excited t 2g levels in 
VOPO to ~ 30 meV. Further experimentation, including infrared spectroscopy and EPR 
measurements, is desirable to locate the excited t2 g orbital levels of V 4+ in VOPO. 

A more likely possibility is that VOPO actually can be modelled accurately as an 
S= 1/2 spin ladder at low energy scales (hence the susceptibility is approximately cor- 
rect), but there are important additional spin interactions that split the one-magnon band 
into two bands, which we observe at 3.5 meV and 6.0 meV at k = n. There are vari- 
ous possibilities for such terms, including ladder-ladder interactions, chain dimerization, 
departures from isotropy, next-nearest-neighbor interactions and so forth. Two examples 
of such interactions are well known. Planar anisotropy in the spin Hamiltonian could 
lead to different gaps for "in-plane" and "out-of-plane" spin-wave modes; this has been 
observed in S= 1 antiferromagnetic chain materials [T1J. A large next-nearest-neighbor 



9 



interaction is believed to be present in the CaV^g spin lattice jT8); similar interactions 
may be present in VOPO as well. 

In conclusion, it appears that a generalization of the simple Heisenberg spin-ladder 
Hamiltonian (1) will be required to describe the magnetic properties of VOPO. The 
correct interpretation of the magnetic excitations would be facilitated by more detailed 
experimental studies of the low-lying bands, in particular by inelastic neutron scattering 
studies of VOPO single crystals. The technical difficulties of growing adequate single 
crystals of VOPO have precluded such studies to date. We hope to carry out such an 
experiment in the near future. 
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Figure Captions 



Figure 1. Predicted spin waves in (VO)2P2C>7 with Jj_ = Jy = 7.0 meV. The one- 
magnon band (solid) and the onset of the two-magnon and three-magnon continua are 
shown. 

Figure 2. Inelastic scattering from (VO) 2 P20 7 at Q — 0.8134 A -1 , at temperatures of 
T=10.8K and 40.8K, measured at HB1A with E< = 14.728 meV and collimation 40'-20'- 
20'- 70'. The gap mode at 3.5 meV and the 6.0 meV mode are clearly visible in the 10. 8K 
data. The corrected neutron counts (see text) are normalized to counts observed in 1250 
monitor units; each monitor unit corresponds to about one second of counting time. 

Figure 3. Constant energy scans of scattering intensity versus momentum transfer 
Q at uj = 3.5 meV (upper panel) and u = 6.0 meV (lower panel), normalized to 500 
monitor units. The dashed vertical lines show the values of Q — nQ n corresponding to 
one-dimensional nuclear (n even) and magnetic (n odd) zone centers. Note the strong 
peaking at Q = Q w in both cases. Measurements were done at HB1A with fixed incident 
energy Ej = 14.728 meV and collimation 40'-40'-40'-70'. 

Figure 4. Constant Q scans at Q — 3Q W and Q = 5Q T at T=1.6K(lower panel) and 
the subtracted difference scattering (upper panel). The scattering intensity is normalized 
to 100 monitor units. The line in the lower panel is a guide to the eye through the 5Q W 
data. The solid line in the upper panel is a least-squares fit of the difference to a sum of 
three Gaussians and a constant background, with fitted centers equal to 3.62(13) meV, 
6.46(16) meV and 14.32(22) meV. The measurements were carried out at HB3 using a 
fixed scattered neutron energy E/ = 14.7 meV and collimation 40'-40'-40'-120'. 
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Figure 1. 
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Figure 2. 
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Figure 3. 
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Figure 4. 
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